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ABSTRACT 

A  model  is  constructed  of  a  repairable  item  inventory  system  where 
on-hand  stocks  are  resupplied  both  by  repair  and  procurement.  The  basic 
properties  of  the  model  are  analyzed.  A  Poisson  failure  process  is 
assumed,  and  steady  state  probabilities  are  found  for  operation  of  the 
model  under  <R5Q>  and  <S-1,S>  inventory  policies.  Further,  a  compound 
Poisson  failure  process  under  an  <S-1  ,S>  inventory  policy  is  also 
analyzed.  Measures  of  supply  performance,  such  as  expected  number  of 
backorders  per  unit  time,  are  derived  for  the  system  under  the  various 
policies. 
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C-,     number  of  items  in  repair,  X(t),  at  t=0 

C2     number  of  items  on  order,  Z(t),  at  t=0 

C3     net  inventory,  Y(t),  at  t=0 

D      delay  between  item  failure  and  establishment  of  a  purchase 
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D,(t)   number  of  items  repaired  and  returned  to  stock  in  the  interval 
(0,tj;  a  random  variable. 


D9(t)   number  of  items  ordered  and  received  in  the  interval  (0,t];  a 
random  variable. 
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H(t)    the  number  of  items  on  hand  (in  stock  minus  backorders)  at  any 
time  t;  a  random  variable. 

I      time  required  to  inspect  an  item;  a  random  variable. 

ICP     Inventory  Control  Point 

I(t)    inventory  position,  i.e.,  on  hand  +  in  repair  +  on  order,  at 
any  time  t;  a  random  variable. 

L      mean  procurement  lead  time. 

L      time  between  placement  of  an  order  and  receipt  of  material 
into  stock;  a  random  variable. 

N,(t)   number  of  failures  entering  repair  during  an  interval  of  length  t. 

Np(t)   number  of  failures  which  are  lost  to  the  system  during  an  interval 

of  length  t. 
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N^  (t)  number  of  failures  placed  on  order  during  an  interval  of  length  t. 

p      probability  an  item  which  has  failed  is  repairable 


p,  probability  a  failed  item  is  returned  to  the  repair  depot. 

p«  probability  an  item  which  arrives  at  the  repair  depot  can 
actually  be  repaired. 

q  probability  an  item  which  has  failed  is  lost  to  the  system. 

Q  the  amount  of  each  order. 

r  probability  a  batch  of  failures  is  repairable. 

R  reorder  level . 

R.  mean  repair  time. 

R0  time  required  to  repair  an  item;  a  random  variable. 

<R,Q>  type  of  inventory  policy  whereby  an  order  is  placed  for  Q  items 
whenever  inventory  position  first  reaches  x  ^  R. 

<s,S>  type  of  inventory  policy  whereby  an  order  is  placed  for  S-x 
items  whenever  inventory  position  first  reaches  x  <_  s. 

T,  in-transit  delay  time,  customer  to  inspection;  a  random  variable, 

T?  in-transit  delay  time,  repair  depot  to  stock  depot;  a  random 

U  number  of  failures  per  batch;  a  random  variable. 

U.  number  of  failures  in  the  i   batch;  a  random  variable. 

V(t)  number  of  batches  of  failures  in  an  interval  (s,s+t]. 

V,(t)  number  of  batches  of  repairable  failures  in  an  interval  (s9s+t]. 

V?(t)  number  of  batches  which  are  lost  to  the  system  in  an  interval 
(s,s+t]. 

W(t)  number  of  items  which  fail  in  an  interval  (s,s+t]. 

X(t)  number  of  items  in  repair  at  any  time  t. 

Y(t)  net  inventory  (on  hand  and  in  repair)  at  any  time  t. 

Z(t)  number  of  items  on  order  at  any  time  t. 

X  rate  of  a  Poisson  process. 
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I.  INTRODUCTION 

Although  inventory  theory  was  one  of  the  first  areas  to  which 
operations  research  analysts  turned  their  attention,  most  of  the  work 
has  been  concerned  with  the  control  of  consumable  items,  i.e.  items 
that  are  consumed  in  use.  Until  fairly  recently  the  control  of 
repairable  items  did  not  excite  a  great  deal  of  effort.  However,  it 
has  now  been  realized  that,  at  least  in  military  applications,  repair- 
able or  recoverable  items  are  a  major  concern  of  inventory  managers, 
both  materially  and  financially.  Repairables  often  account  for  a 
major  portion  of  the  total  dollar  investment  even  though  the  number  of 
items  may  be  small  compared  to  the  number  of  consumable  items.  Further, 
because  repairable  items  are  generally  complex  assemblies,  they  often 
require  special  storage,  handling,  packaging,  and  shipment  methods. 
/  One  of  the  major  difficulties  of  applying  inventory  theory  to 
consumable  items  is  the  extreme  variability  of  item  parameters,  such  as 
demand,  procurement  lead  time,  sources  of  supply,  etc.  Further,  the 
very   fact  of  very   large  numbers  of  items  complicates  an  already  diffi- 
cult task.  Repairable  items,  on  the  other  hand,  have  had  a  tendency 
to  minimize  some  of  these  problems.  The  total  number  of  items  is  much 
smaller  and  it  seems  somewhat  easier  to  determine  accurate  parameters. 
For  instance,  a  large  number  of  repairable  items  may  be  classified 
electronic  sub-components.  Often,  these  can  be  further  grouped  and 
then  treated  as  a  group.  The  sources  of  supply  and  the  sources  of  repair 
are  generally  completely  known  and  limited  in  number.  Repair  lead  times 
and  procurement  lead  times  may  be  often  determined  quite  accurately  with 
limited  variances.  Demand  is  generally  more  stable  in  the  sense  that  if 


a  given  repairable  item  is  a  sub-component  of  a  larger  unit,  it  is 
generally  easy  to  find  out  how  many  larger  units  are  in  use,  and  this 
at  least  puts  an  upper  bound  on  demand. 

On  the  other  hand,  the  control  of  inventories  of  repairable  items 
is  complicated  by  the  fact  of  two  "sources  of  supply",  i.e.,  repair  or 
procurement.  The  purpose  of  this  thesis  is  to  analyze  a  model  of  a 
repairable  item  inventory  system  under  some  common  inventory  policies 
assuming  Poisson  and  Compound  Poisson  demand  processes.  The  question 
of  optimizing  the  system  under  the  various  policies  is  not  taken  up, 
since  inventory  managers  and  other  military  decision-makers  have  not 
yet  agreed  as  to  the  appropriate  optimizing  criteria. 


II.   A  MODEL  OF  A  REPAIRABLE  ITEM  INVENTORY  SYSTEM 

The  model  used  in  this  paper  is  a  fairly  close  approximation  of  the 
actual  repairable  item  inventory  system  in  use  in  the  Navy  today.  Almost 
ewery   item  stocked  by  Naval  vessels  and  installations  is  assigned  for 
technical  and  management  control  to  some  Inventory  Control  Point  (ICP). 
Based  on  economic  and  technical  considerations,  the  ICP  will  designate 
certain  items  as  "repairable."  These  items  are  carried  in  stock,  at 
various  stock  points  or  depots  scattered  throughout  the  Navy  Supply 
System.  Although  for  some  items  the  ICP  does  not  know  how  much  of  an 
item  a  certain  stock  depot  has,  repairable  items  are  almost  always 
"transaction  reporting"  items.  This  means  that  all  issues,  receipts, 
losses,  etc.  of  that  item  at  the  Stock  Depots  are  reported  daily  to  the 
ICP,  wmcn  closely  manages  the  item  at  a  central  level. 

Since  items  are  often  designated  as  repairable  due  to  high  procurement 
cost  or  critically  low  system  stocks,  the  ICP  is  vitally  interested  in 
regaining  all  items  which  have  failed.  To  facilitate  this  return  of 
failed  items,  ships  or  stations  requesting  issue  of  repairable  items  are 
required  to  certify  that  a  like  number  of  failed  items  "has  been  or  will 
be  turned  in."  However,  due  to  the  ship's  operating  position,  distance 
from  stock  points  or  repair  depots,  human  procrastination,  etc.  there  may 
be  some  delay  in  failed  items  arriving  at  the  repair  shop.  Some  items 
may  be  lost  in  transit  or  may  never  be  turned  in.  Those  that  do  arrive 
at  the  repair  depot  are  inspected  and  some  are  found  not-fit-for-repair 
and  are  discarded  for  salvage  value  or  cannabalized  for  usable  sub-parts. 
Those  that  can  be  repaired  are  eventually  returned  to  the  stock  depots 
for  reissue.  Meanwhile  the  ICP,  observing  some  items  being  lost  to  the 


system,  may  decide  to  procure  seme  new  items  to  bring  the  total  system 
stock  back  to  some  desired  level.  The  system  just  described  is  shown 
in  Figure  1 . 
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Figure  1.  Actual  Repairable  Item  Inventory  System 

It  should  be  noted  that  in  actuality  there  usually  will  be  several  stock 
depots  and  probably  at  least  two  repair  depots. 

However,  this  system  can  be  considerably  simplified  with  little  or  no 
loss  of  generality  if  it  is  considered  from  the  point  of  view  of  the  ICP. 
Typically,  an  ICP  will  know  only  when  items  are  issued,  when  items  are 
received  into  stock  from  repair  or  procurement,  and  when  failures  are 
turned  in  to  the  repair  depots.  From  this  information,  the  total  number 
in  system  stock  or  in  repair  at  any  time  can  be  calculated  easily.  The 
ICP  controls  the  procurement  of  new  items  and  thus  knows  the  number  on 
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order  at  any  time.  From  these  remarks,  it  can  be  seen  that  the  various 
delays  in  the  system,  as  seen  from  the  ICP,  are  merely  parts  of  either 
delay  in  repair  or  delay  while  on  order.  Therefore,  it  will  be  assumed 
that  items  which  are  repairable  undergo  a  "repair  time"  which  is  the  time 
from  issue  of  the  replacement  item  out  of  stock  to  receipt  of  the  corre- 
sponding repaired  item  back  into  stock.  Repair  time  in  Figure  1  is 
therefore  T,  +  I  +  R  +  T?.  Likewise,  "procurement  lead  time"  is  the 
time  between  issue  and  return  to  stock  of  a  replacement  for  items  which 
are  lost  or  not  repairable  and  thus  must  be  purchased. 

The  basic  quality  of  failed  items  is  whether  they  are  repairable  or 
must  be  replaced.  Assume  failed  items  are  turned  into  the  repair  depot 
with  probability  p-,  and  of  those  an  item  is  repairable  with  probability 
p?,  then  failed  items  end  up  in  repair  with  probability  p,p?  and  end  up 


i 


'  I  '  / 


Now,  to  simplify  the  model,  assume  that  T-. ,  T^,  I,  and  D  are  equal  to 
zero,  or  equivalcntly  that  the  ICP  has  complete,  instantaneous  knowledge 
of  the  location  and  condition  of  all  items,  and  can  develop  repair  times 
and  procurement  lead  times  as  given  above. 

Let  mean  repair  time  be  denoted  by  R..  Mean  procurement  lead  time 
will  be  denoted  by  L.  Assume  that  R.  has  a  probability  density  function 
fR  and  probability  distribution  function  FR.  Likewise  L  has  density 
f.  and  distribution  F.  .  Let  p,p2  =  p  and  1  -  p,p«  =  q.  Then  the  system 
described  above  can  now  be  modeled  as  shown  in  Figure  2.  Assume  that 
failure  of  any  item  is  independent  of  the  failure  or  non-failure  of  other 
items.  Further  assume  that  the  number  of  customer  requisitions  in  any 
period  of  time  is  a  Poisson  random  variable  with  mean  rate  X.     The  number 
of  items  requested  per  requisition  may  be  one  each,  in  which  case  the 
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Figure  2.  Model  of  the  Repairable  Item  Inventory  System 


process  is  a  simple  Poisson  process.  If  more  than  one  item  can  be 


requestedj  the  procrz",   is  a  compound  Poisson  process 


items  in  repair  at  any  time  t  is  X(t).  The  number  of  items  on  order  at 
any  time  t  is  Z(t).  The  number  on  hand,  which  is  the  number  in  stock 
minus  the  number  of  backorders,  at  any  time  t  is  denoted  H(t).  The 
net  inventory  of  the  system,  defined  as  the  number  in  stock  plus  in 
repair  minus  backorders  at  any  time  t  is  Y(t)  =  X(t)  +  H(t). 
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III.  ANALYSIS  OF  THE  MODEL  UNDER  POISSON  DEMAND 

A.  GENERAL  ANALYSIS  OF  THE  MODEL 

The  number  of  failures  in  any  interval  (s,s+t]  is  assumed  to  be 
Poisson  distributed  with  mean  At.  From  the  assumptions  in  Section  II, 
it  is  seen  that  every   failure  results  in  either  an  additional  item  under- 
going repair  or  an  additional  item  on  order.  Assuming  that  the 
classification  of  failures  are  mutually  independent  and  independent  of 
the  process,  it  has  been  shown  [Parzen  1962]  that  the  number  of  items 
entering  repair  in  the  interval  (s,s+t]  is  Poisson  with  mean  rate  Xp. 
Likewise,  the  number  of  failures  which  will  be  lost  to  the  system  in  the 
interval  (s,s+t]  is  Poisson  with  mean  rate  Xq.  Further,  the  conditional 
distribution  of  the  number  of  items  entering  repair  in  an  interval 

^i ,  i '  i  j    y  i  yen     u::cju    11    iai  'uica    ijc.v^    OCCUrcu     in     Kl\u\*     iiilcl  V"-      10    UiiiOiniu  I 

with  parameters  n  and  p. 

Now  it  can  be  shown  that  the  number  of  failures  entering  repair  in  a 
given  interval  is  independent  of  the  number  of  failures  which  are  lost 
to  the  system.  Let  N,(t)  be  the  number  of  failures  entering  repair 
during  an  interval  of  length  t  and  N2(t)  be  the  number  of  failures  which 
are  lost  to  the  system  during  the  same  interval.  Then  the  joint  proba- 
bility mass  function  for  both  types  of  failures  is: 

PCN^t)  =  x,  N2(t)  "-  n-x) 

■  P(N](t)  =  x,  N(t)  =  n) 

=  PlN^t)  =  x|N(t)  =  n)  P(N(t)  =  n) 


=  (")PX< 
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ni_  px{1.p)n-x  (Xt)x(Xt)n-xe-Pxte-(1-P'U 


x!(n-x)!  n! 

_  px(At)VPAt   (Np)n-x(Xt)n-xe-(1-P)At 
x!  (n-x)! 

.  (PAt)xe'PAt  (gXt)n"xe-qAt 
x!        (n-x)! 

=  P(N1(t)=x)  P(N2(t)=n-x) 

Thus  N,(t)  and  N«(t)  are  independent. 

At  this  point,  before  further  analysis  can  be  done,  the  processes 
of  repair  and  procurement  must  be  given  some  sort  of  formal  definition. 
The  number  of  items  in  the  repair  depot  can  be  viewed  as  the  number  of 
busy  servers  in  a  queue,  specifically  an  infinite  server  queue.  For 
■jt^sus  bp.in"  ropairsd  the  idea  of  an  infinite  serwe^  ou^ue  is  not 
completely  satisfactory.  Certainly  any  repair  depot  set  up  to  handle  a 
continuing  stream  of  failures  will  probably  have  the  ability  to  work  on 
more  than  one  item  at  a  time.  However,  depending  on  the  Poisson  failure 
process  and  the  repair  time  distribution,  the  probability  of  large 
numbers  of  items  entering  or  being  in  repair  at  any  time  may  not  be 
insignificant.  Nevertheless,  the  assumption  will  be  made  that  the 
repair  process  can  be  considered  an  infinite  server  queue,  i.e.  the 
repair  shop  immediately  starts  work  on  all  items  as  soon  as  they  are 
received. 

The  other  specification  necessary  for  analysis  is  the  form  of  the 
reorder  policy.  The  inventory  systems  considered  in  this  paper  are  all 
continuous  review  systems,  that  is,  all  transactions  are  recorded  imme- 
diately, and  thus,  the  inventory  manager  knows  the  exact  value  of  stock 
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on  hand,  on  order,  or  in  repair  at  any  time.  A  common  continuous  review 
policy  is  the  <R,Q>  policy,  which  specifies  that  a  constant  amount  Q  be 
ordered  whenever  the  inventory  position  reaches  R.  The  inventory  posi- 
tion is  defined  as  on  hand  plus  on  order  plus  in  repair  minus  backorders. 
Under  this  policy  the  procurement  process  may  be  thought  of  as  a  queue 
where  the  number  of  busy  servers  at  a  time  t  represents  the  number  of 
orders  of  size  Q  that  are  outstanding  at  that  time.  Another  frequently 
used  policy  is  the  <s,S>  policy,  meaning  order  up  to  S  whenever  the 
inventory  position  falls  below  s.  In  the  case  where  s  =  S-l ,  an  order 
is  placed  everytime  a  failure  cannot  be  repaired.  Under  the  <s,S>  policy, 
the  procurement  process  is  also  comparable  to  an  infinite  server  queue 
where  the  number  of  busy  servers  at  a  time  t  is  the  number  of  orders 
outstanding  at  that  time.  When  items  are  demanded  one  at  a  time,  the 
<[^>o^>  ^clicw  is  the  s^rr!^  ^^  an  ^.S-3,  nolicv. 

Now  let  X(t)  be  the  number  of  items  undergoing  repair  at  time  t  and 
Z(t)  be  the  number  of  items  on  order  at  time  t.  Then 
X(t)  =  ^(t)  +  C1  -  D^t) 

Z(t)  =  Np'Ct)  +  C2  -  D2(t) 

where  C,  is  X(0),  C?  is  Z(0),  D, (t)  is  the  number  of  items  repaired  and 
returned  to  stock  in  the  interval  (0,t],  D«(t)  is  the  number  of  items 
ordered  and  received  in  the  interval  (0,t],  and  N2  (t)  is  the  number  of 
items  lost  to  the  system  which  have  been  placed  on  order  by  the  time  t. 
N2  (t)  is  seen  to  be  related  to  N?(t)  by: 


N2  (t) 


N2(t) 


Q  il 
where  |[nj~  is  the  largest  integer  less  than  n.  If  items  are  reordered 

ewery   time  one  is  lost  to  the  system,  i.e.  Q=l ,  then  N2  (t)  =  N2(t). 
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Otherwise  N2  (t)  <_  N«(t).  Now  both  X(t)  and  Z(t)  can  be  likened  to  the 
number  of  busy  servers  in  a  queue.  D,(t)  depends  only  on  N,(t),  the 
repair  queue  discipline,  and  repair  time  distribution.  N«  (t)  and  D?(t) 
depend  only  on  N2(t),  the  order  queue  discipline  and  procurement  lead 
time  distribution.  Assuming  independence  between  the  two  queue  opera- 
tions, it  is  seen  that  X(t)  and  Z(t)  are  independent. 

Let  H(t)  be  the  number  of  items  in  stock  ready  for  issue  at  time  t, 
where  H(t)  is  negative  if  backorders  exist.  The  net  inventory,  Y(t)  = 
X(t)  +  H(t),  is  the  number  of  items  in  stock  plus  in  repair  minus  back- 
orders.  It  is  important  for  further  analysis,  that  the  independence  of 
X(t)  and  Y(t)  be  shown,  since  H(t)  will  be  obtained  from  X(t)  and  Y(t). 
This  independence  has  been  shown  by  Richards  [Ref.  5]  and  it  is  instruc- 
tive to  repeat  the  proof  here: 

"...exdnrine  the  Pinner  in  which  transitions  are  made  in  the  two 
stochastic  processes,  (X(t),  t  >  0)  and  (Y(t),  t  >_0). 

Whenever  a  type-1  failure  [into  repair]  occurs,  H(t)  decreases  by 
one,  but  X(t)  also  increases  by  one  for  a  net  change  of  zero  in  Y(t). 
On  the  other  hand,  H(t)  decreases  by  one  while  X(t)  remains  constant  when 
a  type-2  failure  [procurement]  occurs;  thus  a  type-2  failure  causes  a  net 
decrease  of  one  in  Y(t)  but  does  not  affect  X(t). 

The  arrival  of  a  procurement  increases  Y(t)  by  the  amount  Q,  while 
the  repair  level  is  unaffected.  Conversely,  X(t)  decreases  while  H(t) 
increases  by  the  like  amount  whenever  an  item  leaves  repair.  Thus, 
repairs  decrease  X(t)  but  leave  Y(t)  unchanged." 
Then  it  can  be  seen  that 

X(t)  =  C1  +  N-,(t)  -  D^t) 

Y(t)  =  C3  -  N'(t)  +  D2(t) 
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where  C3  =  Y(0).  As  previously  shown  D,(t)  depends  only  on  N,(t)  and 

i 
not  on  N2(t).  Conversely,  N2  (t)  and  D2(t)  are  functions  of  N«(t)  but 

not  N-j(t).  Thus  X(t)  is  a  function  only  of  N,(t)  while  Y(t)  is  a  function 

only  of  Np(t).  Since  N,(t)  and  N„(t)  have  been  shown  to  be  independent, 

then  X(t)  and  Y(t)  are  independent. 

B.  AN  <R,Q>  POLICY  WITH  CONSTANT  PROCUREMENT  LEAD  TIME 

Consider  first  the  operation  of  the  model  under  an  <R,Q>  policy, 
where  the  inventory  position  I(t)  =  X(t)  +  H(t)  +  Z(t)  =  Y(t)  +  Z(t). 
An  order  is  placed  for  Q  units  whenever  I(t)  falls  to  R.  Assume  that 
all  delay  times  other  than  repair  time  and  procurement  lead  time  are 
zero.  Further  assume  that  procurement  lead  time  is  a  constant,  L,  Then 
the  "net  demand"  placed  on  the  system  stock,  Y(t),  i.e.  the  number  of 
items  which  must  be  replaced  through  procurement  is  a  Poisson  process  of 
rate  q.  Thus  the  model  from  this  point  of  view  is  merely  an  inventory 
system  operating  under  an  <R,Q>  policy  experiencing  Poisson  demands  at 
rate  q.  This  model  has  been  solved  by  Hadley  and  Whi tin  [Ref.  2]. 
Translating  their  notation  into  that  used  in  this  paper,  the  following 
are  given: 

P[Y(t)=y]  =  1  [F(R+l-y;XqL)  -  F(R+Q+l-y;AqL)]      -~  <y  <  R+l 

=  i  £1  -  F(R+Q+l-y;AqL)]       R+l  <  y  <  R+Q 

oo      i  -r 
where  F(a;r)  =  1  -   E     '•  ■ — 

i=a 

i.e,  the  complementary  cumulative  of  a  Poisson  distribution  of  rate  r. 
Now 

Y(t)  ■  X(t)  +  H(t) 
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thus 

H(t)  =  Y(t)  -  X(t) 

P[H(t)=h]  =  P[Y(t)-X(t)=h] 

oo 

=  I     P[Y(t)=i+h,X(t)=i] 
i=0 

Since  it  has  been  shown  that  Y(t)  and  X(t)  are  independent 

oo 

P[H(t)=hJ  =  E   P£Y(t)=i+h]  P[X(t)=i] 
i=0 

Note  than  an  order  is  placed  for  Q  units  whenever  the  inventory  position 
falls  to  R.  Thus  I(t)  varies  between  R  +  1  and  R  +  Q.  But 

R  +  1  5  I(t)  <  R  +  Q 
implies 

R  -r  1  5  Z(t)  +  Y(t)  <  R  +  Q 

Z(t),  the  number  of  items  on  order,  is  always  non-negative,  i.e.  Z(t)  ^L  0 
Thus 

Z(t)  +  Y(t)  <  R  +  Q 
implies 

Y(t)  £  R  +  Q  -  Z(t)  5  R  +  Q 

or     Max  Y(t)  =  R  +  Q  when  Z(t)  =  0 

Therefore   P£Y(t)=y]  =  0  for  y  >  R  +  Q 

R+Q-h 

means   P[Htt)=h]  =   Z   P[Y(t)=i+h]  P[X(t)=i] 

1=0 

The  probability  that  Y(t)  =  y  has  already  been  developed.  Now  P[X(t)=x] 
is  merely  an  application  of  Palm's  queuing  theorem  which  states  that  for 
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an  infinite  server  queue,  with  Poisson  arrivals  of  rate  *  and  service 
times  distributed  F(t)  with  mean  M,  the  steady  state  probabilities  of  x 
busy  servers  are  Poisson  with  rate  XM,  i.e.  in  this  case, 


PlX(t)=x]  = 


tApRA)VAPRA 


x: 


Since  mean  repair  time  is  taken  to  be  R.,  and  of  course 


E[X(t)]  =  XpRA 


Then 


R+Q-h 
P[H(t)=h]  =   I   P[Y(t)=i+h]  P[X(t)=i] 
i=0 

R+Q-h 


x  '■  l  _ 

2    77  UCR+l-i-h;XqL)"F(R+Q+l-i-h;XqL)]  p(i;XpRfl)   for  i+h  <  R+l 


i=0 
R+Q-h 


where   p(a;r) 


£    ttD  -  F(R+Q+l-i-h;XqL)]  p(i;XpR.)    for  R+l  <  i  +  h  <_  R+Q 
i=0  A 

(r)Vr 


a: 


Simplifying 


R+Q-h 


P£H(t)=h]  =  ~      S  £1  -  F(R-h-i;XqL)  -  1  +  F(R+Q-h-i ;XqL)]  •  p(i;XpRA) 
4  i=0 

l  R+Q"h 
=  7T   ^  [F(R+Q-h-i;Xql_)  p(i;XpR.)  -  F(R-h-i  ;XqL)  p(i;XpRj] 

Q  i=0  A 


r 


-V 


R+Q-h  R+h 

I       p(i;XpRfl)  F(R+Q-h-i;XqL)  -  Z     p(i  ;XpRA)  F(R-h-i;Xql) 
i=0      A  i-0 


Since  H(t)  =  Y(t)  -  X(t)  it  is  easy  to  find  the  expected  value  of  stock 
on  hand: 


19 


E[H(t)j  =  E[Y(t)J  -  E[X(t)J 

Richards  showed  that  the  expected  value  of  net  inventory  for  an  <R,Q> 
policy  is  R  +  ^_ —  m  where  m  is  the  mean  number  of  failures  lost  to 
the  system  in  a  procurement  lead  time.  Thus 

E[H(t)]  =  R  +  S±l  -  XqL  -  XpRA 

=  R  +  311  -  x(pRA  +  qL) 

Now  inventory  position  fluctuates  between  R+l  and  R+Q,  i.e., 

R+l  <  Z(t)  +  Y(t)  <  R+Q 

Assume  the  system  starts  with  Y(0)  =  R  +  Q  and  thus  Z(0)  =  0.  Then  when 
Y(t)  drops  to  R,  Z(t)  becomes  Q.  If  this  order  does  not  come  in  by  the 
time  Y(t)  =  R  -  Q,  then  Z(t)  becomes  2Q  and  so  forth.  ! he  tabic  below 
shows  this  relationship: 

Y(t)  Z(t) 

R  +  Q  0 

*  • 

•  • 

R+l  0 

R  Q 

•  • 

R  -*Q  +  1  Q 

R  -  Q  2Q 

«  • 

R  -  2Q  +  1         2Q 

*  • 

*  * 

•  —  • 

R  -  (n-l)Q         nQ 
R-(nQ)  +  1         nQ 
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Then  Z(t)  =  nQ  if  and  only  if  R-nQ+1  <  Y(t)  <_  R-(n-l  )Q    n=l,2,... 

P[Z(t)=nQ]  =  P[R-nQ+l  <  Y(t)  <  R-nQ+Q] 

=  FYttj(R-nQ+Q)  -  Fy(t)(R-nQ) 

x 
where   F„m(x)  =  £   P[Y(t)  =  v] 

Note  that  if  Q=l ,  i.e.  reorder  every   time  an  item  is  lost  to  the  system, 
this  reduces  to 

P[Z(t)=nJ  =  FY(t)(R-n+l)  -  FY(t)(R-n) 
=  P£Y(t)=R-n+l] 
=  ^  JJ(n;AqL)  -  F(n+l;XqL)J 
iyni  \n„-XqL 

i  A  Vi  L-  i     c 


n! 

That  is,  when  Q=l ,  the  number  on  order  at  any  time  t,  like  the  number  in 
repair,  is  Poisson  distributed  as  expected  by  Palm's  Theorem. 
Now    I(t)  =  Y(t)  +  Z(t) 

E[I(t)]  -  E[Y(t)]  +  E[Z(t)] 

r  +  9+1=  r  +  9+1  .  xqL  +  E[Z(t)] 

or     E[Z(t)>  AqL 

Expressions  have  now  been  developed  for  X(t),  Y(t),  Z(t)  and  H(t)  and 
the  system  is  thus  completely  described.  One  measure  of  effectiveness 
that  is  almost  always  calculated  for  inventory  models  is  the  expected 
number  of  backorders  at  an  arbitrary  point  in  time.  It  is,  of  course 
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00 

E   b  p[H(t)=-b] 
b=l 

where  b  =  the  number  of  backorders.  Under  previous  definition  the 
quantity  on  hand  is  considered  a  negative  value  when  backorders  exist, 
thus  H(t)  =  -b  denotes  b  backorders  on  the  records. 

Another  well-known  measure  of  effectiveness  is  the  expected  number 
of  backorders  per  unit  time  which  is  also  given  by  the  expected  number 
of  failures  minus  the  expected  number  of  fills.  The  expected  number 
of  failures  is  just  A,  the  rate  of  the  Poisson  failure  process.  Now 
a  fill  can  occur  only  if  H(t),  stock  on  hand,  is  greater  than  zero. 
H(t)  >  0  =>Y(t)  -  X(t)  >  0 
=>Y(t)  >  X(t) 

Now  Y(t)  is  never  greater  than  R  +  Q.  If  Y(t)  >  X(t),  but  X(t)  >  0 

(.always),  then  Y(t)  >  0-  So 

p[fillj  =  p[Y(t)  >  X(t)]  for  0  <  Y(t)  <  R+Q 

Therefore: 

R+Q 
p[fillj  =  E  p[Y(t)  >  X(t)|Y(t)=y]  p[Y(t)=y] 
y=l 

p[Y(t)  >  X(t)|Y(t)=y]  =  p[X(t)  <  y] 

y-1  (XPRA)Xe-XpRA 

=  E 

x=0      x! 

=  F(y-l;XpRA) 

Then 

R+Q 
pEfill]  =  E  F(y-1;APRA)  p[Y(t)=y] 

y=i 

R  R+Q 

=  E   p[Y(t)=y]  F(y-l;XpR.)  +   E   p[Y(t)=y]  F(y-1 ;ApRA) 
y=l  y=R+l 
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R  1 

=  Z  -Q-[F(R+Q-y;AqL)  -  F(R-y;AqL)]  F(y-l;XpRA) 

R+Q 
+   E   £  F(R+Q-y;Xql_)  F(y-1  ;ApRj 
y=R+l  g  A 

,  f  r+q  R 

=  ±\     E  F(R+Q-y;AqL)  F(y-l;XpRA)  -  E  F(R-y;.\qL)  F(y-1  ;ApRA) 


Then  the  expected  number  of  backorders  is  X  -  X  •  p[fill] 
R+Q  R 


=  A 


-  q-\  E  F(R+Q-y;AqL)  F(y-1  ;ApRA)  -  E  F(R-y;AqL)  F(y-1  ;ApRA) 


C.  AN  <S-1,S>  POLICY  WITH  STOCHASTIC  REPAIR  AND  PROCUREMENT  LEAD  TIMES 

The  <S-1  ,S>  policy  has  often  been  suggested  as  especially  appropriate 
for  repairable  items  since  the  object  of  the  policy  is  to  maintain  the 
total  system  stock  at  fr,c   level  S.  i.e.  when  an  item  fails  ^nd  cannot  he 
repaired  another  is  immediately  placed  on  order.  Assume  that  there  is 
no  lag  between  a  non-repairable  failure  and  the  placing  of  an  order. 
Then  for  this  model 

X(t)  +  H(t)  +  Z(t)  =  S  for  all  t  >  0 

The  number  in  repair  at  any  time  t,  X(t),  and  the  number  on  order  at  any 

time  t,  Z(t)  are  modeled  as  infinite  server  queues,  with  arrival  rates 

Ap  and  Aq  respectively,  where  q  =  1-p.  In  this  model  both  repair  time 

and  procurement  lead  times  are  random  variables  with  means  RA  and  L 

respectively.  Thus.  Palm's  queuing  theorem  is  applicable  to  both  X(t) 

and  Z(t)  giving 

(ApRfl)VApRA 

P[X(t)=x]  = x=0,l,2,... 

x! 
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P[Z(t)=z]  =  iMjV^!l       z=0,l,2,... 

z: 

i.e.  both  X(t)  and  Z(t)  are  Poisson  distributed. 
Now  note  that  since  X(t)  +  H(t)  +  Z(t)  =  S 

Y(t)  +  Z(t)  =  S 

then    p[Y(t)=y]  =  p[Z(t)=S-y]         y  an  integer,  -»  <_  y  <_  S 

_  (AqL)$-VXqL 

(s-y): 

Also    p[H(t)=h]  =  p[Y(t)-X(t)=h]       h  an  integer,  —  <  y  <  S 


S 
=  I  p[Y(t)=y,X(t)=y-h] 
y=h 

c 

--     E  p[Y(t)=y]  p[X(t)=y-h] 

y=h 


S   (AqL)S^e-X^L  (ApRj^V^A 


p[H(t)  =  h]  =  E 7c — rr-    — ^r 

H       y=h     (s-y):      (y-h): 

Feeney  and  Sherbrooke  [Ref.  4],  using  a  slightly  different  model, 
with  an  <S-1 ,S>  policy,  suggested  the  following  three  measures  of  supply 
performance  (redefined  for  the  present  model): 

(1)  R(S),  called  the  ready  rate,  is  the  steady  state  probability 
that  S  or  less  units  are  in  resupply.  For  the  present  model  "in  resupply" 
means  in  repair  plus  on  order. 
„   R(S)  =  p[X(t)+Z(t)  <  S] 

s 

=  £  p[X(t)+Z(t)=w]    since  both  X(t)  and  Z(t)  are  non-negative 
w=0 
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s 

=  z 

w=0 


w 


Z  p(X(t)=x,Z(t)=w-x) 
x=0 


S  w 

-I  Z  p[X(t)=x]  p[Z(t)=w-x] 

w=0  x=0 

S  w  (ApRfl)xe-APRA   (AqU^e^1 

R(S)  =  Z  Z 

w=0  x=0      xl  (w-x).' 


(2)  F(S),  called  fills,  is  the  expected  number  of  demands  per  time 
period  that  can  be  filled  immediately  from  stock  on  hand.  As  in  the 
last  section  this  is  just  the  failure  rate  times  the  probability  of  a 
fill,  i.e.  X  p[fill] 

p[fill]  =  p[H(t)  >  0] 

=  p[X(t)  <  Y(t)] 


=  £  p[X(t)  <  y|Y(t)=y]  p[Y(t)=y] 

y=i 


=  £  p[X(t)  <  y]  p[Y(t)=y] 
y=l 


S-y  -AqL 


I     F(y-l;APRA)^kLZe 

y-l       A     (s-y)! 


F(S) 


=  \   Z  F(y-l;ApRJ  ■*-*-* 

y=l       A     (s-y): 


Nov/,  it  is  well  known  that  for  a  Poisson  demand  distribution  the  <R,Q> 
and  <s,S>  policies  are  the  same,  with  R  +  Q  =  S  and  Q  =  S-s.  In  this 
case,  with  an  <S-1  ,S>  policy  R  +  Q  =  S,  R  =  S  -  1 ,  and  Q  =  1 .  In  the 
previous  section  it  was  shown  that  for  the  <R,Q>  model, 
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i  \  R+Q  R 

p[fill]  ::  g- \    S  F(R+Q-y;AqL)  F(y-1  ;ApRA)  -  Z  F(R-y;Aql_)  F(y-l;XpRA) 


Substituting  the  values  for  R+Q,  R,  and  Q: 

S  S-l 

p[fill]  =  Z  F(S-y;Aql_)  F(y-l;ApR.)  -  Z  F(S-y-l  ;AqL)  F(y-1 ;ApR  ) 
y=l  A    y=l  A 

S-l 
=  Z  [F(S-y)  -  F(S-y-l);AqL]  F(y-l;ApR,)  +  F(0;AqL)  F(y-l;ApR,) 

y=l  A  A 

s-l 

=  Z  p(S-y;AqL)  F(y-1  ;ApRj  +  F(0;AqL)  F(y-1  ;ApR J 
y=l  A  A 


Now  F(0;AqL)  =  p(0;AqL)  =  p(S-S;AqL). 

Thus 

S 
p[fill]  -  Z  p(S-y;AqL)  F(y-l;ApRj 
y=l  A 

which  agrees  with  the  expression  given  above  for  the  <S-1,S>  policy. 
The  expected  number  of  backorders  in  a  unit  time  is,  again, 

A  -  A  •  p[fill] 

S 

=  A  -  A  Z  F(y-l;ApRA)  p(S-y;AqL) 

y=i      A 

(3)  A  third  measure  of  supply  performance  considered  is  S(S),  the 
number  of  "units  in  service,"  defined  as  the  expected  number  of  units 
in  routine  resupply  at  a  random  point  in  time.  Routine  resupply  is 
further  defined  as  the  number  of  units  in  resupply  minus  the  number  of 
units  backordered.  Again  "in  resupply"  is  taken  to  mean  in  repair  plus 
on  order.  We  note  that  if  X(t)  +  Z(t)  =  r  5  S,  then  H(t)  _>  0,  i.e., 
there  are  no  backorders,  and  the  number  of  units  in  routine  resupply  is  r. 


26 


On  the  other  hand  if  X(t)  +  Z(t)  =  r  >  S,  HCt)  <  0,  i.e.,  there  are 

backorders,  and  the  number  of  units  in  routine  resupply  is  S. 

Thus 

S 
S(S)  =  Z  r  P[X(t)+Z(t)=r]  +  Z  S  P[X(t)+Z(t)=r] 
r=0  r=S 

S   f  r  ")        °°    f  r 

=  Z  r-Jz  P[X(t)=w]  p[Z(t)=r-w]l  +  Z   S  4  Z  P[X(t)=w]  p[Z(t)=r-w] 
r=0   )w=0  (   r=S    )w=0 


S         (   r     (APRfl)VXpRA     (AqL)r-VAqL 

S(S)  =     Z     r\     Z 

r=0       )  w=0  w!  (r-w)l 


-     (  r      (ApRjVApRA   (XqL)r-We-^L 

+  s      z    J    z ■ 

r=S  )w=0  wl  (r-w)! 
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IV.  ANALYSIS  OF  THE  MODEL  UNDER  COMPOUND  POISSON  DEMAND 

The  model  in  the  last  section  assumed  that  failures  were  a  Poisson 
process  of  rate  A,  i.e.  the  number  of  failures  in  a  time  interval  of 
length  T  has  a  Poisson  probability  mass  function  with  parameter  AT. 
This  implied  failures  occurred  one  at  a  time,  and  thus  each  requisition 
for  replacement  was  for  exactly  one  item. 

In  practice,  the  Poisson  model  may  be  quite  reasonable  for  major 
repairable  units,  i.e.  items  of  high  cost  and/or  large,  complex  structure. 
However,  there  are  smaller  not  so  expensive  items  that  are  also  economi- 
cally repairable  and  that  may  experience  failures  in  batches,  i.e.  more 
than  one  at  a  time.  One  good  example  might  be  the  printed  circuit  "cards" 
found  in  many  modern  computers,  radars,  and  other  electronic  equipment. 
These  cards  are  generally  easily  removed  and  installed,  and  under  condi- 
tions of  peak  load,  current  surges,  etc.,  several  may  fail  at  the  same 
time.  Generally,  the  only  damage  may  be  to  the  resistors,  capacitors, 
transistors,  etc.  which  are  soldered  onto  the  printed  circuit,  while  the 
circuit  card  itself  is  undamaged  and  repair  is  feasible. 

These  types  of  items  lead  to  the  consideration  of  failure  processes 
in  which  the  number  of  failures  in  a  given  time  period  is  a  compound 
Poisson  process  (W(t),  t>0}.  This  process  can  also  be  represented  by 

V(t) 
W(t)  ■  2   un 
n=l 

where  (V(t),  t>0}  is  a  Poisson  process  which  is  the  number  of  customer 
requisitions  received  and  {Un,  n=l,2,...}  is  a  family  of  independent 
identically  distributed  random  variables  where  Un  is  the  number  of  items 
requested  or  the  nth  requisition.  It  is  necessary  that  (V(t),  t>0}  and 
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U   are  assumed  to  be  independent  and  for  this  model  that  means  that 
the  number  of  items  per  requisition  is  independent  of  which  requisition 
it  is,  which  is  an  entirely  reasonable  assumption.  It  is  noted  that  the 
simple  Poisson  process  is  just  a  special  case  of  a  compound  Poisson 
process,  where  the  probability  {U  =l,for  any  n}  =  1 .  The  compound 
Poisson  is  an  extremely  rich  family  of  distributions.  Various  seemingly 
quite  different  processes  can  be  described  merely  by  the  proper  choice 
of  a  compounding  distribution. 

Assume  that  W(t)  the  number  of  items  which  fail  in  a  unit  time 
interval  is  a  compound  Poisson  process  with  rate  A  and  that  U  is 
distributed  as  a  random  variable  U,  with  probability  mass  function  f 
and  characteristic  function  <J>.  .  Then  it  has  been  shown  [Parzen,  p.  130] 
that  the  characteristic  function  of  U(t)  is 

a    ( Y\  -  pxtvfi|'vx;~i ; 

Mtr  '  - 

Referring  to  the  repairable  item  inventory  model  of  this  paper, 

assume  that  all  failure  items  in  a  particular  batch  are  either  repairable 

or  all  are  non-repairable  and  result  in  a  procurement.  Further,  let  the 

probability  that  all  items  in  a  batch  are  repairable  be  r.  Then  the 

number  of  items  to  enter  repair  from  the  n   batch,  R  ,  is 

U  with  probability  r 

0  with  probability  1-r 
That  is    prob  [R  =  u]  =  prob  [U  =  u]  r     u=l,2,3,... 

prob  [R  =  0]  =  1-r 


The  characteristic  function  of  R  is 


*R  =  EreieRn]  .  i    eieo  p[R a 
j=o 


-  P[Rn=o]  ♦   i    ei*j  p[       3 


00 

=  1-r  +  r     Z     e1<})j  p[U  =j] 
j=l  n 

=  1-r  +  r4>y 

n 


^R      =  "I   -  r  +  r^u 


Then  N, (t),  the  number  of  items  entering  repair  in  a  time  interval    (0,t] 

is 

V(t) 
N,(t)  =     I       R 
1  n=l       n 

and  from  above 


i . .  >   1  \ 


\(t)ix)  =e        Rn     ;"' 
=  gAtd-r+r^-l) 

=  eXtr(4>u-l) 

which  is  just  the  same  compound  Poisson  process  as  failures  with  rate  X 
replaced  by  rate  Xr.  Likewise,  N?(t),  the  number  of  items  which  are 
placed  on  order  in  an  interval  (0,t],  can  be  shown  to  be  the  same 
compound  Poisson  process  with  rate  X(l-r). 

Now,  let  the  number  of  batches  entering  repair  be  V-,(t)  and  the 
number  of  batches  which  must  be  replaced  by  procurements  be  V?(t),  so 
that  V(t)  =  V-(t)  +  Vp(t).  The  number  of  batches  in  the  compound 
Poisson  model  corresponds  to  the  number  of  items  in  the  simple  Poisson 
model.  In  Section  II  it  was  shown  that  the  number  of  items  entering 
repair,  N, (t),  and  the  number  of  items  lost  to  the  system,  N2(t)s  were 


independent  random  variables.  In  the  compound  Poisson  process,  this 
implies  that  V-,(t)  and  V2(t)  are  independent.  Let  N^t)  be  the  number 
of  failures  in  an  interval  (s,  s  +  t)  which  can  be  repaired,  and  let 
N  (t)  be  the  number  of  failures  in  the  same  interval  which  are  lost  to 
the  system.  Then 

V(t) 
N  (t)  =  E    UnXl 
n=l 

V(t) 
N  (t)  =  Z        Unx2 
n=l 


_f|  if  the  n 

cl  JO  otherwis 

fo  if 

x2  JJ  otherwise 


th 


batch  goes  to  repair 
where      Xn  Mn  nthprwise 

the  nth  batch  goes  to  repair 


By  assumption  [Un,  n=1,25...}  is  a  Family  oi  indcj.enflei 3 , 

distributed  random  variables,  and  a  batch  is  either  repairable  or  lost 
to  the  system.  Thus  N^t)  and  N2(t)  are  sums  of  independent  identically 
distributed  random  variables,  and,  furthermore,  no  Un  which  is  in  the 
sum  for  N^t)  is  also  in  the  sum  for  N2(t)  and  vice  versa.  Therefore, 
N  (t)  and  N2(t)  are  independent.  Further,  X(t),  the  number  of  items  in 
repair  at  any  time  t  is  a  function  only  of  N^t)  and  the  repair  queue 
parameters.  Likewise,  Z(t),  the  number  of  items  on  order  at  any  time  t 
is  a  function  only  of  N2(t),  the  procurement  policy,  and  the  procurement 
queue  parameters.  Thus  X(t)  and  Z(t)  are  still  independent  for  the 
compound  Poisson  case,  assuming  a  batch  goes  completely  to  repair  or  is 
completely  lost  to  the  system. 

Again  defining  H(t)  to  be  the  number  of  items  on  hand  at  time  t,  and 
calling  Y(t)  ■  X(t)  +  H(t)  the  net  inventory  at  time  t,  the  same  analysis 
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as  in  Section  II  shows  that  X(t)  and  Y(t)  are  independent.  Using  the 
notation  of  that  section 

X(t)  =  C1  +  I^Ct)  -  D^t) 

Y(t)  =  C3  -  N2'ct)  +  D2(t) 

where  C-,  =  X(0),  C3  =  Y(0),  and  N~  (t)  is  the  number  of  items  which  have 
been  ordered  in  the  interval  (0,t].  It  has  been  shown  that  N,(t)  and  fL(t) 
are  independent  and  it  is  assumed  that  D-,(t)  is  a  function  only  of  N,(t), 
i.e.,  that  the  number  of  items  repaired  in  an  interval  of  length  t  is  a 
function  only  of  the  number  of  items  which  enter  the  repair  queue  and  the 
queue  operation  itself,  and  that  N«  (t)  and  D?(t)  are  functions  of  N?(t). 
Thus  X(t)  and  Y(t)  are  functions  of  independent  random  variables  and  are 
themselves  independent  random  variables. 

Assume  that  the  inventory  ^y~tr~,   Is  governed  by  an  <S-1  ,S>  policy 
whereby  every   failure  that  cannot  be  repaired  is  immediately  placed  on 
order.  This,  of  course,  may  result  in  an  order  for  more  than  one  item 
under  the  assumption  that  an  entire  batch  is  either  repairable  or  not. 
Assume  further  that  all  items  which  enter  repair  or  which  are  placed  on 
order  at  the  same  time  have  the  same  repair  time  or  procurement  lead 
time  respectively.  The  effect  of  these  assumptions  is  that  items  which 
fail  in  the  same  batch  are  either  repairable  and  are  repaired  as  a  batch, 
or  that  those  which  are  not  repairable  must  be  ordered  and  are  received 
into  stock  as  a  batch.  While  these  assumptions  may  seem  unduly  restric- 
tive, in  actual  practice  they  are  quite  reasonable.  Certainly  items 
ordered  at  the  same  time  are  generally  received  at  the  same  time.  Further, 
repair  often  consists  of  a  series  of  operations  to  be  performed.  It  is 
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not  unusual  that  the  first  operation  is  completed  on  all  items  in  a 
batch  before  the  second  operation  is  started,  etc.,  and  upon  completion 
of  repair  the  items  are  returned  to  stock  as  a  batch. 

Let  R.  be  mean  repair  time  and  L  be  the  mean  procurement  lead  time. 
Let  f  (x)  be  the  "compounding"  probability  mass  function,  that  is 
fu(x)  =  P{U  =  x}  . 

Now  Feeney  and  Sherbrooke  have  shown  that  under  the  above  assumptions 
given  a  compound  Poisson  arrival  density  of  rate  a  into  a  queue,  the 
steady  state  probability  of  x  units  in  the  queue  is  the  same  compound 
Poisson  density  with  rate  XT,  where  T  is  the  mean  service  time  in  the 
queue.  Since  the  rate  of  items  entering  repair  is  Xr: 

Prob  {x  units  in  repair  at  any  time} 

=  Prob{X(t)  =  x} 


\  ■  •  •  ■  V"  '   - 


=  E fu  U) 

i=0      i!      u 

i* 

where  f   is  the  i-fold  convolution  of  f  . 

Likewise,  with  the  rate  "into"  procurement  being  X(l-r), 
Prob{z  units  on  order  at  any  time} 
=  Prob{Z(t)  =  z} 

-  I     [A(l-r)L]VX(1-r>L  fi*(z) 
i=0        i! 

Then  assuming  the  <S-1  ,S>  policy,  where  an  item  is  ordered  every   time  it 
is  lost  to  the  system,  the  inventory  position,  which  is  equal  to  on  hand 
+  in  repair  +  on  order  -  backorders,  is  always  equal  to  S.  In  terms  of 
the  model 

H(t)  +  X(t)  +  Z(t)  =  S 
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Then 


or  Y(t)  +  Z(t)  =  S 

Prob  (Y(t)=y) 
=  Prob(S-Z(t)=y} 
=  Prob(Z(.t)  =  S-y} 


S  y  [Ad-DLl^^^  fi*  t   , 


z 

1=0 


i : 


Furthermore,  since  Y(t)  =  H(t)  +  X(t) 


H(t)  =  Y(t)  -  X(t) 


-°°  <  H(t)  <  S 


Then    Prob  (H(t)=h} 

=  Prob(Y(t)-X(t)=h> 

S-h 
=  I  Prob(Y(t)=y+h,X(t)=y} 
y=0 

S-h 
=  Z     Prob{Y(t)=y+h>  Prob(X(t)=y> 
y=0 

and  since  Y(t)  and  X(t)  are  independent 


S-h 

Z 

y=0 


S-y-h     [AO-rKlV^1-^ 


Z 
1=0 


i: 


y  (XrRA)VXrRA     ]ie 

f  (S-y-h)  Z__^ f^  (y) 

u       i=0 


The  system  is  now  described  for  any  time  t  >  0.  Again,  the  measures  of 
supply  performance  suggested  by  Feeney  and  Sherbrooke  are  ready  rate, 
fills  and  units  in  service. 

Ready  rate,  R(S)  is  the  steady  state  probability  that  S  or  less  units 
are  in  repair  or  on  order,  i.e.  X(t)  +  Z(t)  <_  S. 

R(S)  =  Prob  [X(t)  +  Z(t)  <  S] 


=  E  Prob  [X(t)  +  Z(t)  =  w] 
w=0 


S  w 

=  Z  E   Prob  [X(t)=x,Z(t)=w-x] 

w=0  x=0 

S  w 

I  Z   Prob  [X(t)=x]  Prob  [Z(t)=w-x] 

w=0  x=0 


S       w 

x 

Z       Z 

I 

w=0  x=0 

i=0 

x  (xrRfl)V*rRA 


^A 


W 


i: 


u    i=0 


-x  [XO.r)L]1e-X(l-r,L 


fjVx) 


Fills,  F(S),  is  the  expected  number  of  demands  per  time  period  for  an 
item  that  can  be  filled  immediately  from  stock  on  hand.  Feeney  and 
Sherbrooke  reasoned  that,  "There  will  always  be  one  fill  if  S-l  or   less 
units  are  in  resupply  [in  repair  plus  on  order]  since  each  customer 
places  an  order  for  at  least  one  unit.  There  will  be  a  second  fill  if 
S-2  or  less  units  are  in  resupply  and  the  customer  order  is  for  two  or 


"II.'!  J-1   'i^i.ui  •  .«....  .  . 


Therefore,  F(S)  is  the  expected  number  of  customers  per  time  period, 
i.e.  X,  times  the  expected  number  of  fills  per  customer.  As  shown  above, 
this  can  be  represented  as 

F(S)  =  A  ^  (probability  a  customer's  demand  is  greater  than  0)  ■ 

[probability  X(t)  +  Z(t)  <_  S-l) 

+   (probability  a  customer's  demand  is  greater  than  1)  • 

(probability  X(t)  +  Z(t)  <  S-2) 

+ 

» « >, « • 

+       [probability  a  customer's  damand  is  greater  than  S-l)    • 
(probability  X(t)  +  Z(t)  <   0.)  V 

=    x[l    •    R(S-l)   +  [l-fu(l)]R(S-2)  +   [l-fuU)-   fu(2)]R(S-3)   + 
,..  +  [l-fuO)-fu(2)-...-fu(S-l)j   R(0)"j 
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This  is  the  same  formula  given  by  Feeney  and  Sherbrooke,  except  the  ready 
rate  has  been  adjusted  for  the  double  "sources"  of  supply,  i.e.  repair 
and  procurement.  Then  the  expected  number  of  backorders  in  a  unit  time 
is,  with  probability  one,  the  expected  number  of  demands  less  the  expected 
number  of  fills,  i.e. 

E[number  of  backorders]  =  XE[U]  -  F(S) 

The  other  measure  of  supply  performance  is  units  in  service,  S(S),  which 
is  the  expected  number  of  units  in  routine  resupply  at  a  random  point  in 
time.  There  are  r  items  in  routine  resupply  if  X(t)  +  Z(t)  =  r  <^  S,  and 
there  are  S  items  in  routine  resupply  if  X(t)  +  Z(t)  >  S.  Thus 


S(S)  =  I     r  Prob[X(t)+Z(t)=r]  +   E   S  Prob[X(t)+Z(t)=r] 
r=l  r=S+l 


l     r  E   Prob£X(t)=w]  Prob[Z(t)=r-w] 

r=l   w=0 


+  Z         S   E   Prob[X(t)=w]  Prob[Z(t)=r-w] 
r=S+l    w=0 


S     r 

E  r   E 

r=l    w=0 


'w    (ArRA)VXrRA       r-w  [X(l-r)L]1e"x"-r'L  „ 

Z    ^ f]   (w)   E —f"    (r-v 

i=0      i!     u     i=0       i:      u 


oo    r 
+  S   E     E 
r=S+l  w=0 


w   (ArRA)VXrRA  .* 

—  fu  (w) 


i=0  i! 

r-w  [x(l-r)L]ie-X(1'r)L   .* 

E  —  f  ,  (r-w) 

i=0        i! 
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V.  CONCLUSIONS 

The  results  developed  in  this  paper  are  relatively  straightforward 
extensions  of  models  and  results  found  by  Hadley  and  Whitten  and  by 
Feeney  and  Sherbrooke.  The  model  is  a  fairly  good  approximation  of 
reality.  The  assumption  of  Poisson  or  compound  Poisson  demand  processes 
has  found  much  favor  in  inventory  theory.  The  steady  state  probability 
of  the  number  of  items  in  repair,  in  stock,  or  on  order,  have  been 
found,  under  various  inventory  policies.  From  these  probabilities, 
formulas  for  several  measures  of  supply  performance  have  been  derived. 
Other  measures  can  be  developed  in  a  similar  manner.  Thus,  given  the 
steady  state  system  probabilities  and  the  derived  measures  of  effective- 
ness, inventory  managers  can  now  attempt  to  optimize  the  system  under 
various  criteria. 

Further  analysis  to  improve  the  model  may  also  be  possible;  for 
instance,  the  probability  of  a  fill  given  a  delay  time  t  in  which  to 
make  the  fill  would  be  a  useful  result.  It  may  also  be  possible  to 
relax  the  assumption  that  under  compound  Poisson  demand,  failure  batches 
either  go  entirely  to  repair,  or  are  lost  to  the  system.  It  would  seem 
more  natural  that  each  item  of  each  batch  may  go  into  repair  or  may  be 
lost  to  the  system. 
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